Ballistic spin currents in mesoscopic metal/In(Ga) As/metal junctions 
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We investigate the ballistic spin transport through a two-dimensional mesoscopic 
metal/semiconductor/metal double junctions in the presence of spin-orbit interactions. It is 
shown that real longitudinal and/or transverse spin currents can flow in the presence of the Rashba 
and Dresselhaus terms. 
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I. INTRODUCTION 

Since the advent of "spintronics" to utilize electron's 
spin rather than its charge for information processing and 
storage,^ there has been growing interest in generating 
spin currents in diverse ways^i^*^i^ Though injecting 
spin-polarized carriers electrically remains a challenge;^ 
various kinds of all-semiconductor devices using ferro- 
magnetic semiconductor heterostructureS^ or spin-orbit 
(SO) interactions^ have been proposed. The so-called 
extrinsic spin Hall effect due to SO dependent scatter- 
ing from magnetic impurities manifests a spin current 
because of a transverse spin imbalance generated from 
a charge current circulating in a paramagnetic metal.'* 
In the weak impurity scattering limit but with substan- 
tial SO couplings, on the other hand, it was suggested 
that the intrinsic spin Hall effect gives rise to a dissipa- 
tionless spin current perpendicular to the external elec- 
tric fieldi^ Moreover, the spin Hall conductance has a 
universal value. However, it was demonstrated^ that the 
dissipationless (unreal) spin current does not vanish even 
in thermodynamic equilibrium in the absence of external 
fields, putting the interpretation of intrinsic spin Hall ef- 
fect in a controversy. 

In this paper we study the ballistic spin transport 
through a mesoscopic double-junction system consist- 
ing of a semiconductor stripe sandwiched by two normal 
metal leads; see Fig. ^ We use the coherent scattering 
theory and show that in the presence of SO couplings, 
both longitudinal and transverse spin currents can flow 



in the semiconductor, 
are real; see Ref. 



It is stressed that these currents 
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II. MODEL AND SCATTERING THEORY 

We consider a two-dimensional electron system (2DES) 
of semiconductor (S) between two normal (N) metal 
leads. We choose such a coordinate system that x-axis 
(y-axis) is perpendicular (parallel) to the N/S interfaces 
and 0-axis is perpendicular to the 2D plane; see Fig. ^ 
The length (width) of the semiconductor is L (W); we 
will consider the limit W ^ oo. Within the effective- 
mass approximation^i the Hamiltonian reads as 



2 m[x) 



v + v{x,y) + nR{x) + nD{x). (1) 



FIG. 1: A schematic of the system. 



The position-dependent effective mass m{x) has values 
of TTie and m* = e,„me in the normal metals and the 
semiconductor {—L/2<x<L/2), respectively. The con- 
finement potential has a potential barrier of height Vq 
inside the semiconductor: 

Vix, v) = [e(.T + L/2) - e{x - L/2)] + Viy), (2) 

where 0{x) is the Heaviside step function and V{y) ac- 
counts for the finite width W. The potential barrier 
height Vq is lower than the Fermi energy Ep in the nor- 
mal metals so that Ep = Ep — Vq > 0. The RashbaiS^ and 
Dresselhaus^- SO coupling terms are given by 

T-Lr = ^{a-xPy-a-yPx) and Hd = ^(o-yPy-a^Px): (3) 

respectively, inside the semiconductors while they vanish 
in the normal metal sides. In Eq. (Q, cr = {cr^, Cy, Cz) 
are the Pauli matrices. 

The Rashba term Tin arises when the confining poten- 
tial of the quantum well lacks the inversion symmetry, 
while the Dresselhaus term Tiu is due to the bulk inver- 
sion asymmetry. In some semiconductor heterostructures 
(e.g., InAs quantum wells) Hr dominatesii, and in oth- 
ers (e.g., GaAs quantum wells) Ti^) is comparable to (or 
even dominant over) Tin^^ . The coupling constants may 
range around a ~ O.leV-A and (3 ^ 0.09 eV-A, respec- 
tively, depending on the structure and material^ 
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FIG. 2: Fermi contours (solid/dotted lines for jj, = ±), and 
on each of them, the wave vectors k (black solid arrows), the 
group velocities (dark/light green arrows), and the spin 
orientations ifik (dark/light blue arrows) of the eigenstates. 
(a) a = 0.5(3 and (b) a = /3 with Ep = 14meV, /? = 0.1 eV-A, 
and Em — 0.063. 



Inside the semiconductor, electrons feel a fictitious, in- 
plane magnetic field in the direction 



nk = X cos ifk + y sm ipj^ . 



(4) 



where t/^k — arg[(/3A;2: — aky) + i{akx — f^ky)]. Accord- 
ingly, the eigenstates with spin parallel (/x = and 
anti-parallel {fi = —) to Ak for a given wave vector 
k — fc(xcos0 -|- ysini/)) are written in the spinor form 



g-l-iVk/2 
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The corresponding energies are 

ft2 



2m* 



[k' ~ 2^lkso{4>)k] 



(5) 



(6) 



where kso{(t)) = (jn%lf?)\J c? ^ — 2Q;/3sin20. From 
the continuity equation for the charge density, one can 
get the expression for the charge current associated with 
a given wave function ^'(r)^'^ 

je = eRe[*^(r)v*(r)] , (7) 

where v is the velocity operator defined by 

^ ^ 'm*~li ^^^^ ~ '^''^^^ ~ li ^^^^ ~ '^'^^^ ■ 

In the same manner, we define the spin currenlii 

■ i-\ ^Tt/ N v(n • cr) + (n • cr)v 
J.(n) = 2* 

according to the continuity equation 

dtQ. + V-\s = Ss (10) 

for the spin density (with respect to the spin direction n) 



g,(n)^-[vl/t(r)(n.^)*(r)] (11) 



and the spin source 



5.(n)--Re 



*^(r)- [7^,n-cr]*(r) 



(12) 



The appearance of the spin source term in Eq. (|10|l is 
not surprising because the spin-orbit couplings break spin 
conservation inside the semiconductor. 

Before going further, it will be useful to understand the 
origin of the spin current in physical terms. As illustrated 
in Fig. 12 for a, /? 7^ the Fermi contours, 

fc^^(0) = /xfc,„(0) + y^fc2^(0) + fcj? (13) 



with k*p = 2m* E*p j are no longer isotropic, and 

the group velocities = ^'j^V'J'j^ of the eigenstates in 
Eq. (|3J) are not parallel to the wave vector Nev- 
ertheless, Eq. (|13|l reveals an important symmetry prop- 
erty of the group velocities: vr|^+ = v~ . It means that 

the two eigenstates with opposite spin orientations make 
the same contributions to the charge transport along the 
k direction (and opposite contributions along the the per- 
pendicular direction) . On the contrary, for the spin trans- 
port with n = rij^, the inverse occurs: two eigenstates 
make opposite contributions along the k direction and 
same ones in the perpendicular direction. This implies 
the possibility of the observation of the net spin current 
flowing perpendicular to the charge current. Particularly 
interesting are the cases of a = ±/3, where all the spin 
orientations, ifik, for different wave vectors are parallel 
or anti-parallel to each other (iy9k = t^I^'i see Fig. |2(b). 
It results from conservation of [ux^Oy] /V2, and the spin 
state become independent of the wave vectoriiiife 

Now we study charge and spin transport in N/S/N 
double junction structures. Coherent scattering formal- 
ism at the N/S interfaces has been thoroughly developed 
in the previous studiesr^*i^ considering the Rashba SO 
effect and appropriate boundary conditions requiring the 
conservation of probability current normal to the inter- 
face. It is straightforward to extend the scattering theory 
to incorporate the Dresselhaus effect. We have used the 
transfer-matrix formalism to calculate the conductance 
through jind inside the semiconductor: for details refer 
to Kefs. ilTiiTl 

We will consider electrons incident from the left lead 
and reflected from the junction interfaces or transmitted 
through them to the right lead. The wave vector of the 
incident electron is at angle 6 with respect to the nor- 
mal to the interface; see Fig. ^ Contrary to the Rashba 
effect, the Dresselhaus effect is not invariant under the 
rotation, leading to anisotropic transport»i2. Hence the 
relative orientation, ^, of the crystal symmetry axes and 
the interface (see Fig. ^ affects especially the spin cur- 
rent signiflcantly. Below we will calculate the charge con- 
ductance G'^u\e) ( v — x,y) in the I'-direction for a defl- 
nite incident angle 9 as well as the angle-averaged quan- 



tity Gi'^ 



d%d9 Gl'^'ie). Also calculated are the 



Ac) 
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(analogously defined) spin conductances Gi^''^\o) and 
Gi'*'"'' polarized in the direction n. The typical values 
for the parameters we will use below are Ep — 4.2 eV, 
e„^ = 0.063, P = O.leV-A, L = 200 nm, and W = 1 fim. 
a ranges from —2/3 to +2f3, and Ep from to 20meV. 
We assume sufficiently low temperatures (fcsT <C Ep). 



III. NORMAL INCIDENCE 



Owing to the symmetry 



[see the dis- 



cussion below Eq. H12() ]. for normal incidence {0 — 0} the 

charge current is completely longitudinal; i.e, G^\6 = 
0) = 0. For a single transverse mode, we obtain the 
longitudinal charge conductance 



32k2 



h 1(1- 



-(1-K)2e2*^'=i|2 



(14) 



where Afc = V^so(-0 + ^ 



*2 
F ' 



Ak/ejrikp, kp = 



y/2meEp/h. Moreover, the spin current has only trans- 
verse component and is polarized entirely in the xy-plane; 



i.e., Gi^'")(0 = O) 
The ns 
by 



for any n and Gy''^\ 



-polarized spin conductance Gy 



' = 0) = 0. 
= 0) is given 



' = 0) = 



(l+«2) 



'kpkso{-C) 

.2\ sin 2AfeL 



(1- 



-K 



2AkL 



|(l-hK)2 - (l-K)2e2»AfeL|2 



(15) 



G 



0) and g|,''"*^(6I = 0) are plotted in Fig. El as 



functions of Ep and a/ (3 for different crystal orientations 

^. The peaks in Gi''\e = 0) and g1,"'"*^(6' = 0) as a 
function of Ep come from the Fabry-Perot interference, 
which gives rise to resonances for 



AkL = ' 



(" = 0,1, 2,. 



(16) 



Unlike the (longitudinal) charge current, the spin current 
is very sensitive to the SO coupling strengths, a and /3, 
and the crystal orientation, ^, as seen from the factor 
Q!;5cos^ in Eq. ifT^ . 



IV. ANGLE-AVERAGED CONDUCTANCES 
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FIG. 3: The charge conductance Gi''\e^O) [(a) and (b)] and 

the spin conductance Gy"'"*' (6 = 0) [(c) and (d)] for normal 
incidence as functions of Ep [(a) and (c)] and a//3 [(b) and 
(d)]. a/ 13 = 0.5 in (a) and (b), and E*p = 12 and 14meV in (b) 
and (d). ^ has been chosen to be (black line), jg (red solid 
line), J (green dashed line), and j (blue dotted line). Notice 
that in (a) curves for different ^'s overlap almost completely. 
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FIG. 4: Angle-averaged charge conductance G^f ' as a function 
of (a) Ep and (b) a//3. ^ = has been chosen, but Gi'^' is 
not sensitive to 5. Values of other parameters are indicated 
in the figures. 



Apparently, the longitudinal charge current has a main 
contribution from the normal incidence. Consequently, 
as shown in Fig. 0] the 0-averaged longitudinal conduc- 

(c) . 

tance Gx is rather similar to the normal incidence case 
gI'^'' (9 = 0), although the peaks are rounded off. 

This is not the case for the spin transport. Figure [S] 
shows the ^-dependence of the spin conductances polar- 
ized in the Ax and z, respectively. Again, the peaks cor- 



For true one-dimensional (ID) leads {kpW ^ 1), 
where only a single transverse mode is allowed, one has 
only to consider normal incidence (9 = 0); or at a certain 
fixed 9^^. In the opposite limit (i.e., kpW 00), there 
are many transverse modes contributing to the transport. 
In this case, we should add up all the contributions from 
6 in the range (— 7r/2, 7r/2). It is quite complicated (even 
though possible) to find the scattering states for non- 
zero incidence angle 9, and more convenient to work nu- 
merically. Therefore, here we just present the numerical 
results. 




FIG. 5; Angle dependences of the spin conductances 
G$,"'"*'(6I) and G$,"'*'(e) for E*p = 14meV, a/P = 0.5, and 
5 = 0. 
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FIG. 6: Transverse spin conductance gJ,^'^' as a function of 
(a) E*p for a/P = 0.5 and (b) a/f3 for £J = 14meV. Inset: 
Gy"'^'' as a function of /3/a with a = O.leV-A fixed. The 
values of ^ are the same as in Fig. |3] 

respond to the Fabry-Perot-type resonances. When sum- 
ming up, the contributions to the fix-polarized spin cur- 
rent from different angles are mostly canceled with each 
other, and hence the angle-averaged spin conductance 
Gy'*'""^ becomes very small compared with the longitu- 

(c) 

dinal charge conductance Gx . On the other hand, the 
z-polarized spin current is not subject to such cancella- 
tions, and remains rather large (still smaller than the lon- 
gitudinal charge current); see Fig. [HI Especially for ^ = 0, 

the spin conductance Gy^'^^ remains almost constant in 
the region \a/(3\ < 1 and changes its sign abruptly at 
a/ (3 = ±1. This behavior is reminiscent of the intrinsic 
spin Hall conductances in the previous worksi^ How- 
ever, in our case Gy*'^^ depends on the strength of the 
SO couplings, the potential barrier, crystal orientation, 
and the channel length, showing no universal character- 
istics. 

Here we stress the differences between origins in the 
spin Hall conductance of ours and the intrinsic spin Hall 
effect. The intrinsic spin Hall effect is an (semi-classical) 
effect driven by external electric field penetrating the (in- 
finitely large) systemi^ In our case, the external bias 
voltage merely shifts the relative chemical potentials of 
the "contacts" (or reservoirs) attached to the metallic 
leads where the electrons undergo ballistic transport and 
does not feel an electric field »iS Moreover, it has been 
pointed out that the spin current in the intrinsic spin 
Hall effect is an equilibrium background current and is 
not real.^ To the contrary, our spin currents originate 



from non-equilibrium properties of the system, and are 
real. 

We also note that when electrons are incident oblique 
to the junction interface, the transverse charge current 
and the longitudinal spin current do not vanish any 

(c) 

longer. Therefore, the angle-averaged conductances Gy 
and gI'''"'' are finite, even though quite small compared 

(c) (c) 

to Gx . Finite Gy in the semiconductor can be at- 
tributed to the anisotropy introduced by the Dresselhaus 
effect. It distorts the group velocity, which thus prefers 
one of ±y directions so that the current has same sign as 
af3. Nonzero Gx^'"'' reflects the breaking of spin conser- 
vation inside the semiconductor. For oblique incidence, 
there exists no direction consistent with the boundary 
conditions along which spin state is stationary (for ex- 
ample, the spin parallel to the direction nt for a given 
wave vector k = kxX + kyy is not stationary any longer 
for k' = —k'^x + ky-y after reflection from the junction 
interface). This means that an electron with any spin 
polarization experiences precession during transmission 
through the semiconductor. 



V. CONCLUSION 

Ballistic spin currents with different spin polarizations 
through mesoscopic metal/2DES/metal double junctions 
have been investigated in the presence of spin-orbit inter- 
actions. Using the coherent scattering theory we showed 
that longitudinal and/or transverse spin currents can flow 
through 2DES. It was argued that arising from the non- 
equilibrium distribution of electrons, the spin Hall cur- 
rents observed are real. 



Acknowledgments 

M.L. thanks W. Belzig, C. Bruder, and J. Schlicmann 
for helpful discussions. This work has been supported 
from the SKORE-A program and the eSSC at Postech. 
M.-S.C. acknowledges the support from KIAS, where 
part of the work was done. 



1 G.A. Prinz, Science 282, 1660 (1998); S.A. Wolf, D.D. 

Awschalom, R.A. Buhrman, J.M. Daughton, S. von 

Molnar, M.L. Roukes, A.Y. Chtchelkanova, and D.M. 

Treger, Science 294, 1488 (2001). 
^ P.R. Hammar, B.R. Bennett, M.J. Yang, and M. Johnson, 

PRL 83, 203 (1999). 
^ Y. Ohno, D.K. Young, B. Beschoten, F. Matsukura, H. 

Ohno, D.D. Awschalom, Nature 402, 790 (1999). 

M.I. Dyakonov and V.I. Perel, Zh. Eksp. Ter. Fiz. 13, 657 

(1971) [JETP 33, 467 (1971)]; J.E. Hirsch, Phys. Rev. 

Lett. 83, 1834 (1999); S. Zhang, Phys. Rev. Lett. 85, 393 



(2000); L. Hu, J. Gao, and S.Q. Shen, Phys. Rev. B 68, 

115302 (2003); 68, 153303 (2003). 
^ S. Murakami, N. Nagaosa, and S.C. Zhang, Science 301, 

1348 (2003); J. Sinova, D. Culcer, Q. Niu, N.A. Sinitsyn, 

T. Jungwirth, and A.H. MacDonald, Phys. Rev. Lett. 92, 

126603 (2004); J. Schliemann and D. Loss, Phys. Rev. B 

69, 165315 (2004). 
® N.A. Sinitsyn, E.M. Hankiewicz, W. Teizer, and J. Sinova, 

cond-mat/0310315 S.-Q. Shen, cond-mat/0310368 
^ E.I. Rashba, Ph ys. Rev. B 68, 241315 (2003); 
|co"nd-mat/0404723"| 



5 



G. Bastard, Wave Mechanics Applied to Semiconductor 
Hetero structures (Halstead, New York, 1988); M.G. Burt, 
Phys. Rev. B 50, 7518 (1994). 
^ Y.A. Bychkov and E.I. Rashba, J. Phys. C 17, 6039 (1984). 

G. Dresselhaus, Phys. Rev. B 100, 580 (1955). 

B. Das, S. Datta, and R. Reifenberger, Phys. Rev. B 41, 
8278 (1990); G. L. Chen, J. Han, T. T. Huang, S. Datta, 
and D. B. Janes, Phys. Rev. B 47, 4084 (1993); J. Luo, 

H. Munekata, F. F. Fang, and P. J. Stiles, Phys. Rev. B 
41, 7685 (1999). 

G. Lommer, F. Malcher, and U. Rossler, Phys. Rev. Lett. 
60, 728 (1988); B. Jusserand, D. Richards, H. Peric, and 
B. Etienne, 69 p. 848 (1992); B. Jusserand, D. Richards, 
G. AUan, C. Priester, and B. Etienne, Phys. Rev. B 51, 
4707 (1995). 

J. Schliemann, J.C. Egues, and D. Loss, Phys. Rev. Lett. 



90, 146801 (2003). 
^'^ L.W. Molenkamp, G. Schmidt, and G.E.W. Bauer, Phys. 
Rev. B 64, 121202 (2001); G. Feve, W.D. Ohver, 
M.Aranzana, and Y. Yamamoto, Phys. Rev. B 66, 155328 

(2002) ; 

T. Matsuyama, C.-M. Hu, D. Grundler, G. Meier, and U. 

Merkt, Phys. Rev. B 65, 155322 (2002). 

J. Schhemann and D. Loss, Phys. Rev. B 68, 165311 

(2003) . 

Gonzalo Usaj and C.A. Balseiro, ljcond-mat/0405065 
V. Marigliano Ramagha, D. Bercioux, V. Cataudella, G. 
De Filippis, and C.A.Perroni, cond-mat/0403534 (unpub- 
lished) . 

S. Datta, Electronic Transport in Mesoscopic Systems 
(Cambridge Univ. Press, Cambridge, 1995). 



